Mon. Not. R. Astron. Soc. 000, [TH5] (0000) Printed 19 January 2012 (MN WF&L style file v2.2) 



One step beyond: The excursion set approach with 
correlated steps 



Marcello Musso 1 ' 2 * k Ravi K. Sheth 1 ' 3 



(N 

o 

(N 
C 

00 



o 

U 
9 L1, 

6 



1 77ie AMits Salam International Center for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy 

2 CPS, Universite de Louvain, 2 Chemin du Cyclotron, 1348 Louvain-la-Neuve, Belgium 

3 Center for Particle Cosmology, University of Pennsylvania, 209 S. 33rd St., Philadelphia, PA 19104, USA 



19 January 2012 



> 

oo 
cn 

o 

(N 



ABSTRACT 

We provide a simple formula that accurately approximates the first crossing distribu- 
tion of barriers having a wide variety of shapes, by random walks with a wide range 
of correlations between steps. Special cases of it are useful for estimating halo abun- 
dances, evolution, and bias, as well as the nonlinear counts in cells distribution. We 
discuss how it can be extended to allow for the dependence of the barrier on quantities 
other than overdensity, to construct an excursion set model for peaks, and to show 
why assembly and scale dependent bias are generic even at the linear level. 

Key words: large-scale structure of Universe 



1 INTRODUCTION 

In hierarchical clustering models, to estimate cluster abun- 
dances at any given time one must estimate the abundance 
of sufficiently overdense regions in the initial conditions 
(Press & Schechter 1974). The problem is to find those re- 
gions in the initial conditions that are sufficiently overdense 
on a given smoothing scale, but not on a larger scale. The 
framework for not double-counting smaller overdense regions 
that are embedded in larger ones is known as the Excursion 
Set approach (Epstein 1983; Bond et al. 1991; Lacey & Cole 
1993; Sheth 1998). 

In this approach, one looks at the overdensity around 
any given random point in space as a function of smoothing 
scale. The resulting curve resembles a random walk, whose 
height tends to zero on very large smoothing scales. In this 
overdensity versus scale plane, the critical density for col- 
lapse defines another curve, which we will call the barrier. 
The double-counting problem is solved by asking for the 
largest smoothing scale on which the walk first crosses the 
barrier. It is fairly straightforward to solve this problem nu- 
merically, by direct Monte-Carlo simulation of the path in- 
tegrals (Bond et al. 1991). This is particularly simple if the 
steps in the walk are independent, but one can also include 
correlations between the steps, whose nature depends on the 
underlying fluctuation field (i.e., for a Gaussian field, on the 
power spectrum) and the form of the smoothing filter. 

When the steps are independent, exact solutions for 
constant (Bond et al. 1991) or linear (Sheth 1998) barriers 
are known. Exact solutions for more general barrier shapes 
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are not known, but good analytic approximations are avail- 
able (Sheth & Tormen 2002; Lam & Sheth 2009). In the 
appropriate units, these solutions are self-similar, i.e. inde- 
pendent of the form of the power spectrum (of course, they 
depend strongly on the barrier shape). However, an exact 
solution of the first crossing problem for correlated steps 
is still unknown. The main goal of the present work is to 
provide a simple formula which works for a wide variety of 
barrier shapes, smoothing filters, and power spectra. This is 
done in Section [2] - our main result is equation JSJ , and it 
is explicitly not self-similar. Section [3] describes a number of 
extensions of this calculation - having to do with walks con- 
ditioned to pass through a certain point in a certain way, or 
with barriers whose height depends on hidden variables. A 
final section summarizes our results, indicating how we ex- 
pect our work to be used when fitting to the halo abundances 
which, recent simulations indicate, are not quite self-similar. 



2 FIRST CROSSING DISTRIBUTION WITH 
CORRELATED STEPS 

In what follows, we will assume that the underlying fluctu- 
ation field is Gaussian. We comment on non-Gaussian field 
in the Discussion section. In hierarchical models, the vari- 
ance of the fluctuation field is a monotonic function of the 
smoothing scale (the exact relation depending on the shape 
of the power spectrum and the smoothing window.) There- 
fore we can use the terms smoothing scale and variance in- 
terchangeably, and we will use s to denote the variance. 

Let B(s) denote the height of the barrier on 'scale' s, 
and p(8, s) the probability that the walk has height 8 on 
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this scale. We will assume that (<5) = 0, so that s = {S 2 }. We 
would like to write down the probability f(s) that 5 < B(S) 
for all S < s and 5 > B(s) at s. When the steps in the walk 
are uncorrelated the two conditions separate, simplifying the 
analysis. For correlated steps, this simplicity is lost. 



2.1 The completely correlated limit 

Recently Paranjape et al. (2012) have argued that there is 
considerable virtue in thinking of the limiting case in which 
the steps in the walk are completely correlated. In this case, 
if the walk had height So on scale So, then it has height S = 
So \Js/ So on scale s, so 5 is a smooth monotonic function of s. 
Therefore, if 5 first exceeds B(s) on scale s, it was certainly 
below B(s) on all S < s, and one need not account for this 
requirement explicitly. Hence, the first crossing distribution 
is just 



/cc(s) = 



dSp(S, s) 



B(s) 



cxp(-g 2 (s)/2s) d£(s)/\/5 
ds 



for S < S m in, (1) 



where 5* m i n is the scale on which B(s)/^ r s is minimum (in 
many cases of interest Smin = oo). Note that this relates the 
shape of f(s) to that of p(8, s) on the same scale. 

Paranjape et al. (2012) show that, despite the strong as- 
sumption about the deterministic smoothness of the walks, 
this expression provides a very good description of the first 
crossing distribution (at small s) even when the steps are 
not completely correlated. Physically, this is because if one 
thinks of real walks as stochastic zig-zags superimposed on 
smooth completely correlated trajectories, at small s a walk 
has not fluctuated enough to depart significantly from its de- 
terministic counterpart. One might thus expect that there 
is danger of double counting trajectories with two or more 
crossings only when s becomes large. 

In cosmology the large s regime is not nearly as in- 
teresting as the small. The discussion above suggests that 
it will be useful to construct an expansion in terms of the 
number of times walks cross the barrier. The first step in 
this program is to assume that no walks double-cross, but 
the actual correlation structure scatters their crossing scale 
around the completely correlated prediction. Accounting for 
this fact alone should allow to estimate f(s) with a greater 
regime of accuracy than equation {l}. Accounting for one 
earlier crossing should be even more accurate, and so on. 



2.2 A bivariate approximation for the strongly 
correlated regime 

To proceed, we assume that when steps are strongly corre- 
lated one can replace the requirement that S(S) < B(S) for 
all S < s, which is a condition on all the steps in the walk 
prior to s, with the milder requirement that S(s — As) < 
B(s — As) for As — > - a condition on the one preceding 
step. Because now the walk values on the two scales are inde- 
pendent, the analysis is more involved than when the steps 
were completely correlated, but the increase in complexity 
is relatively minor because we only require a bivariate dis- 
tribution. 

For small As we can expand both S and B in a Taylor 



series. The condition on the walk height at the previous step 
means that 5(s) — AsdS/ds < B(s) — As dB /ds. If we use 
primes to denote derivatives with respect to s, then the first 
crossing distribution of interest is given by the fraction of 
walks which have 

B(s) < 5 < B(s) + As (8' — B') and S'>B'(s). 

That is to say, 

rB(s)+As (S'-B 1 ) 



f(s) ds — lim 



dd' 



dSp(6, 5') 



B(s) 



Asp(B,s) dS'p(S'\B)(S' - B') 



(2) 



For a Gaussian field p(B, s) is Gaussian, and the con- 
ditional distribution in the integral is also Gaussian, with 
shifted mean and reduced variance. If we define 



(SS') 2 



the mean of p(6'\B) is (S'\B) = 7 B{ S' 2 }^ 2 /( S 2 ) l ' 2 and 
the variance is ( S' 2 )(1 — j 2 ). (Our notation was set by the 
fact that, for Gaussian smoothing filters, 7 equals the quan- 
tity spectral quantity a 2 /ao(T2 which Bardeen et al. 1986 
called 7. For tophat filters, the integrals over the power spec- 
trum which define our 7 are given by Paranjape et al. 2011.) 
Note that (86') = 1/2, and thus (S'\B) = B/2s. 

Before we evaluate the integral, note that the com- 
pletely correlated approximation corresponds to the limit 
in which 7 = 1 and p(S'\B) becomes a delta function 
centered on B/2s. The integral then yields B/2s — B' = 
— y r sd(B/y r s)/ds, and the resulting expression is consistent 
with equation {TJ. Notice that our analysis has indeed ex- 
tended the completely correlated solution by replacing the 
delta function with a Gaussian whose width depends explic- 
itly on the underlying power spectrum. 

In the generic case, one still has an integral over a single 
Gaussian distribution. If we define 



S'-B' 



P(s) 



B(s) 



and /3,(s 



-m 



din f3(s) 
din v/s 



(the reason for our choice of sign for /3» will become clear 
later), then equation @ gives 



sf(s) 



dx x ■ 



-( 3: -7/3.(s)) 2 /2(l-7 2 ) 



27 V2tt 

Evaluating the integral yields 



sf(s) 



,-/3 2 M/2 

2^/2^ 



■P. 



1 + erf (r^/^2) e 



-r 2 /3j/2 



27rr/3» 



(4) 



(5) 



(recall /?* depends on s); this is our main result. 

Comparison with equation ([T]) shows that the term in 
square brackets above represents the correction to the com- 
pletely correlated solution. There are two points to be made 
here. First, this correction term depends on V, indicating 
that the shape of the first crossing distribution depends ex- 
plicitly on the form of the underlying power spectrum. In 
this respect, walks with correlated steps are fundamentally 
different from walks with uncorrelated steps. Second, when 
/3« 3> 1 then this term tends to unity so the first crossing 
distribution reduces to that for completely correlated walks. 
To see what large /3* implies, it is useful to consider some 
special cases. 
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Figure 1. Distribution of the scale y = 8 2 /s on which walks 
first cross the barrier <5 C (1 + as/8 2 ). Histograms show numerical 
Monte Carlo results for Gaussian smoothing of P(k) oc k~ 1 - 2 , 
for which T 2 = 9/10; top to bottom are for barriers with a = 
— 1,0 and 0.5. For comparison, the two dotted curves show the 
corresponding distributions for a = and —1 when steps are 
uncorrclatcd. Smooth curves show Eq. © with T 2 = 9/10 and 
the appropriate values of a; the agreement indicates that our 
formula works well for a wide variety of barrier shapes. Symbols 
with error bars show results for a constant barrier (a = 0) and 
tophat smoothing of a ACDM P(k). This shows that, in contrast 
to when steps are uncorrelated, the first crossing distribution does 
indeed depend (weakly) on P(k). The solid curve shows Eq. (0 
with 7 = 1/2 for which T 2 = 1/3; the agreement shows that our 
formula works well for a wide range of T. 



2.3 When the barrier has constant height 

For a constant barrier, B(s) = S c and /3*(s) = 8 c /^/s. The 
latter is conventionally denoted as v, so that ft* 3> 1 corre- 
sponds to large v. In terms of v, equation ((5} is 



erf(r^/y / 2) + c 



2-kYv 



(6) 



In this form, it is clear that T~ acts as the i^-scale below 
which the correction to the completely correlated distribu- 
tion becomes important. (Note however that our Y is not 
the same parameter defined by Peacock & Heavens 1990 
and used in Paranjape et al. 2012.) 

Figure Q] compares this formula to distributions gen- 
erated via Monte-Carlo simulation of the walks (following 
Bond et al. 1991). We present results for two rather different 
power spectra and smoothing filters, with two different val- 
ues of P. Our first choice is Gaussian smoothing of P(k) oc k n 
with n = -1.2, for which F 2 = (n+3) /2 = 9/10. Our second 
is tophat smoothing of a ACDM power spectrum; in this case 
7 itself depends on scale, with 7 w 1/2 (and hence F w 1/3) 
on the scales where v ~ 1. The figure shows that equa- 
tion ([6]), with r 2 = 9/10 and 1/3 respectively, provides an 
excellent approximation to the Monte-Carlo distributions. 



Note in particular that, for ACDM, ignoring the scale de- 
pendence of r (by simply using 1/3 on all scales) works very 
well. This demonstrates that Eq. © is a simple and accurate 
approximation for the CDM family of models. 

2.4 Linear barriers 

For linear barriers B(s) = <5 C (1 + as/<5 2 ), making /3(s) = 
v-\-a/v and /3* = v — ol/v. Figure [1] shows results for a — — 1 
and 1/2, both for Gaussian smoothing of P(k) oc fc -1 ' 2 ; 
Equation © with F 2 = 9/10 works very well. Note in par- 
ticular that our formula is quite different from the Inverse 
Gaussian distribution associated with uncorrelated steps. 



3 EXTENSIONS 

The analysis above has been so simple, and its results so 
accurate, that it is interesting and natural to extend it to a 
variety of other problems, some of which we outline below. 



3.1 Halo bias 

The excursion set approach is often used to quantify the 
correlation between halo abundances and their environment. 
This is done by computing the ratio of f(s\8o, So), the first 
crossing distribution subject to the additional constraint 
that the walks passed through some So on some large scale 
5*0 before first crossing S c on scale s ^> So, to f{s). 
Motivated by the previous Section, we now set 



f(s\So,S )ds= / dS' 



in the limit where As 



S(s)+As (S-B 1 ) 

dS P (5,5'\8o,So) (7) 

B(«) 

0. This will give 



f(s\6 , So) = p{B\8 , So) / dS'p(5'\B, So) (5' - B') , 

J B' 



(8) 



which is very similar to equation @, only with modified 
mean values and variances. The Gaussian outside the in- 
tegral has mean (<5|5o) = So{ SSo }/So and variance s(l — 
£ 2 ), where £ = (S5 ) /Vs3o~. The one inside has mean 
(S'\B,S ) = ( Ms B + At <So)/(l - £ 2 )> where Ms = [(SS') - 
(S'8 )(8So)/S]/s and Mo = [(80S') - (S5')(S5 ) /s]/S . The 
first Gaussian obeys the scaling assumed by Paranjape et 
al. (2012), who showed it was a good approximation to their 
conditional Monte Carlo distributions. Therefore, it is inter- 
esting to see if this scaling also holds for the integral. 

If {5q8') ~ (or, more carefully stated, if the term 
containing this quantity is smaller than the other) then 
fi B = {58') 1 8 and no = -{SS'){SS ) /sS , so (8'\B,S ) = 
Hb{B — (<5|(5o)). Since this is the same rescaling of B as 
for the first Gaussian, in this approximation f(s\5a,Sa) fol- 
lows the scaling with So that was assumed by Paranjape et 
al. (2012). Therefore, the conclusions of Paranjape & Sheth 
(2012) about the difference between real and Fourier space 
bias factors also hold for our calculation. In particular, the 
Fourier space bias b\ will be simply given by differentiating 
In sf(s) with respect to B(0), whereas the real space cross 
correlation will carry an additional factor of (SSo) /So- For a 
constant barrier 61 = d\ni> f(v) / d8 c yields 



bi = 



1 



1 + 



1 + erf (Tv/y/2) V^Fv 



-r 2 v 2 /2 



(9) 
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The first term on the right hand side is the bias associ- 
ated with completely correlated walks. The correction term 
is negligible for v ^> 1, and it tends to 1/5 C when 

The expression above assumes that (S'So) is negligible. 
This is reasonable when So >C s; e.g. for Gaussian smooth- 
ing filters, {5' So) = a 2 x /2a 2 falls to zero faster than (S5o). 
But at intermediate scales it is not, and /(s|(5o, So) wm have 
additional dependence on So. In Musso, Paranjape & Sheth 
(2012) we show that this will generically introduce scale de- 
pendence into the Fourier space bias factors, even at the 
linear level. 



3.2 Conditional mass functions 

The excursion set theory can also be used to compute halo 
progenitor mass functions and merger rates (Lacey & Cole 
1993). For this, we need the joint distribution f(s,b,S,B) 
of walks that first cross the barrier B on scale S, and first 
cross the barrier b > B on scale s > S. Dividing by f(S, B), 
which we already have, will give the conditional distribution 
/(s, b\S, B). The most natural way to estimate it is to set 

i poo rB(S)+AS(A'-B') 

f(s,b\S,B)^— dA' / dA 

^ s J B> Jb(S) 

l-oo Ms)+As(6'-b') n (X X' A A'l 



f{S)dS 



P(B) 
f(S) 



dA' (A' - B 1 ) p(A'\B) p{b\B, A') 



x / dS' (5' -b')p(S'\b,B,A') 



P{B) 
f(S) 



dA'(A'-B')p(A'\B)f(s\B,A'). (10) 



(Strictly speaking, we should adjust the limits on the inte- 
grals over A to ensure that A < b.) We have written the 
final expression in a suggestive form, to show that the left 
hand side should be thought of as a weighted average over 
first crossing distributions, each with its own value of A'. 

In the limit where the scales are very different, s ^> S, it 
should be a good approximation to ignore the {SA'), {5' A) 
and (S'A') correlations. This makes b-{b\B, A') = b-£(B- 
GA')/(1 — G 2 ) where G is the same quantity as 7, but on 
scale S. This means there is some dependence on A'. Note 
that if there were no correlation with A' (in effect G — > 0) 
then the conditional distribution would have the same form 
as the unconditional one, except that b — > & — BS X / S, and 
s — > s—S(Sx / S) 2 , where we have defined S x = {SA) . This is 
similar to the rescaling for sharp-fc filtering, for which S x = 
S. However, for Gaussian smoothing of a power law P(k), 
G — 7, so one may not set G — > if one does not also set 
7 — > 0. Therefore, the case with G^O may be more relevant. 
In this case {S'\b, B, A') = ■y(l — G 2 )(b~(b\B, A')). Since this 
same term appears in the exponential of p(b\B, A'), we have 
that f(s\B,A') has the same form as equation ((5j), except 
for the shift b — ?> b — {b\B, A'). Therefore, except for the erf 
piece, the final integral over A' can also be done analytically. 
We leave a comparison of this approximation, and the full 
expression (in which we have not ignored (<5'A) etc.), with 
the numerical Monte-Carloed distributions, to future work. 

That f(s\B, A') depends explicitly on A', even in this 
approximation, is significant, since it shows that the con- 



ditional distribution for first crossing b depends not just on 
the fact that B was crossed, but on how B was crossed. If we 
interpret 'how B was crossed' as a statement about the mass 
distribution on smoothing scales larger than S, then the fact 
that f(s\B, A') depends on A' indicates that the formation 
history of the mass within S depends on the surrounding 
environment. This shows that our formalism will naturally 
give rise to 'Assembly Bias' effects of the sort identified by 
Sheth & Tormen (2004), and studied since by many others. 



3.3 Dependence of B on parameters other than 8 

In triaxial collapse models, the barrier is a function of the 
values of the initial deformation tensor (rather than just its 
trace). This makes 

/pe p p poo 

de J dp J de' J dp J dS' 

r5 c (e,p)+As (S'-B 1 ) 

x / dSp(e,p, S, e',p, S') (11) 

J S c (e,p) 

where B' = e' (dB /de)+p' (dB/dp). If the scale dependence 
of e and p can be neglected, then this will simplify to 

/(*) = J de J dpp(e,PA(e,J>)) 

x / dS'S'p(S'\e,p,S c (e,p)). (12) 
Jo 

Further, if the correlation between 8' and e and p can be 
neglected, then 

f(s) = J dep dpg(e,p\8 c (e,p))f(s\8 c (e,p)), (13) 

and the first crossing distribution is that for S c (e,p), 
weighted by the probability of having e and p. If we use 
equation (A3) of Sheth et al. (2001) for the distribution of e 
and p given S, then our analysis should be thought of as gen- 
eralizing their equation (9). In particular, because the result 
is a weighted sum of first crossing distributions, it exhibits 
exactly the sort of stochasticity discussed in Appendix C of 
Paranjape et al. (2011). Performing this calculation more 
carefully is the subject of ongoing work. 

3.4 Excursion set model of peaks 

One of the great virtues of our approach is that it provides 
a simple way to compare the excursion set description with 
that for peaks. The relation is particularly simple for Gaus- 
sian smoothing filters, since then 7 of our equation @ is 
the same parameter which plays an important role in peaks 
theory. This correspondence means that our parameter x is 
essentially the same as the peak curvature parameter; the 
only difference is that we define the derivative with respect 
to the variance s, whereas peaks theory derivatives are with 
respect to smoothing scale R. This means that our inte- 
grals over S' are really just integrals over curvature, which 
makes intuitive sense. Hence, to implement our prescription 
for peaks, we only need to account for the fact that the dis- 
tribution of curvatures around a peak position differs from 
that around random positions. If we write our equation Q 
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as (/3»/2) e 13 / 2 /y/2n times {x)/"/f3„, then we need only re- 
place (x) — > (xf(x)), where f(x) is given by equation (A15) 
of Bardeen et al. (1986). 

For a constant barrier, our 7/3* equals their x#, so our 
expression is simply their equation (A14) weighted by x/x* 
and integrated over x. (Their additional factor of {2-k) s,2 rI 
is just the usual ra/p factor in excursion set theory, which 
converts from mass fractions to halo abundances.) Omit- 
ting the x/x* factor when performing the integral over x 
yields the usual expression for peaks (their equation A18). 
Therefore, one might think of their (A18) as representing the 
'completely correlated limit' for peaks, whereas our analy- 
sis yields what should be thought of as the moving barrier 
excursion set model for peaks. 



4 DISCUSSION 

We presented a formula, equation J5J, which provides an 
excellent description (see Fig. [1} of the first crossing distri- 
bution of a large variety of barriers, by walks exhibiting a 
large range of correlations (i.e., it is valid for a wide class of 
power spectra and smoothing filters). As we discuss below, 
we expect a special case of it - equation ([6]) - to be a good 
physically motivated fitting formula for halo abundances. 

We then showed that our approach provides a simple 
expression for the first crossing distribution associated with 
walks conditioned to pass through a certain point, and hence 
a simple expression for how halo bias factors are modified 
because of the correction term (equation [5]) . We sketched 
why a generic feature of this approach is that even the lin- 
ear bias factor should be scale dependent. We also showed 
how to approximate the first crossing distributions associ- 
ated with two non-intersecting barriers: the probability that 
a walk first crosses barrier B on some scale S and then b > B 
on scale s > S (equation llOp . This exhibit 'assembly bias' so 
they may provide useful approximations for halo progenitor 
mass functions and merger rates. 

Finally, we argued that our approach makes it par- 
ticularly easy to see how the first crossing distribution is 
modified if the barrier depends on hidden parameters, such 
as those associated with the triaxial collapse model (Sec- 
tion ° r peaks (Section 13. 4|l . In our approach, peaks dif- 
fer from random positions only in that the integral over x in 
equation ((4]) is modified. A similar integral in the expression 
for peak abundances leads to scale dependent bias even at 
the linear level (Desjacques et al. 2010), and is why we now 
expect ^-dependent bias even for random positions. 

We are not the first to have considered the corre- 
lated steps problem. Peacock & Heavens (1990) identified 
T (our equation [3]) as the key parameter. Their approxima- 
tion for /(s) is more accurate than more recent approxima- 
tions (Maggiore & Riotto 2010; Achitouv & Corasaniti 2011) 
which are, in any case, restricted to special combinations of 
power spectra and smoothing filter (Paranjape et al. 2012). 
However, our equation © is simpler and more accurate for 
a wider range of barrier shapes, power spectra and smooth- 
ing filters than any of these previous studies, and it can be 
easily extended. Besides the problems outlined above, an ob- 
vious direction would be to include non-Gaussianity, along 
the lines of Musso & Paranjape (2012). This extension is con- 



ceptually straightforward - equations ((2J will hold also for 
a non-Gaussian field - and is currently being investigated. 

We argued that this accuracy stems from the small s 
behavior of the walks. In this regime, which is of most inter- 
est in cosmology, f(s) is well-described by equation {1} for 
completely correlated walks (Paranjape et al. 2012). These 
walks are deterministic: the distribution of their heights is 
a delta function. For real walks, instead, it has a width that 
increases with s. Equation (© was obtained by allowing for 
this broader distribution of heights in the one step prior 
to the crossing. However, it does not explicitly account for 
walks which may have criss-crossed the barrier more than 
once. Musso et al. (2012) discuss how accounting for more 
zigs and zags may yield even greater accuracy at larger s. 

Equation (0 is explicitly the completely correlated first 
crossing rate times a correction factor. It depends on power 
spectrum only because this factor depends on F. Since this 
factor is small at small s, in this regime (which is the usual 
one in cosmology) one should expect to see departures from 
self-similarity, but they should be small. E.g., for a barrier 
of constant height, the first crossing distribution becomes 
equation ([6]), and the correction factor becomes important 
at v < F _1 , with T" 1 ~ 3 for the ACDM family of P(k). 

In cosmology the first crossing distribution is often used 
as a fitting formula. For this purpose, one might treat either 
S c or r or both as free parameters (even though equation [6] 
describes the first crossing distribution very well with no 
free parameters). In this case, the value of S c will depend on 
a variety of factors (Sheth et al. 2001; Maggiore & Riotto 
2010b; Paranjape et al. 2012), and we expect F to depend 
on the effective slope of the power spectrum. In this sense, 
our formula provides a simple way to understand, interpret 
and quantify the departures from self-similarity which sim- 
ulations are just beginning to show. 
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